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Addressing the feasibilty of quantum communication with electrons we consider entangled spin 
states of electrons in a double-dot which is weakly coupled to in-and outgoing leads. We show that 
the entanglement of two electrons in the double-dot can be detected in mesoscopic transport and 
noise measurements. In the Coulomb blockade and cotunneling regime the singlet and triplet states 
lead to phase-coherent current and noise contributions of opposite signs and to Aharonov-Bohm and 
Berry phase oscillations in response to magnetic fields. These oscillations are a genuine two-particle 
effect and provide a direct measure of non-locality in entangled states. We show that the ratio of 
zero-frequency noise to current (Fano factor) is universal and equal to the electron charge. 
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Entanglement and non-locality of Einstein-Podolsky- 
Rosen (EPR) pairs Q| are remarkable features of quan- 
tum mechanics which give rise to striking phenomena 
such as violation of Bell inequalities and secure quantum 
communication 0|. Such and related phenomena have 
been tested with great success for photons but not 
yet for massive particles such as electrons, particularly 
in a solid state environment. There are two immediate 
problems. First we need a scheme by which entangle- 
ment of electrons can be generated in a controlable man- 
ner. In the context of quantum computing || we have 
recently shown that such a scheme can be realized in 
tunnel-coupled quantum dots each of which contains one 
single (excess) electron whose spin defines the qubit |i]-[| ■ 
There are several motivations for such a qubit scheme, 
most notably very long spin decoherence times and 
the qubit defined as electron-spin is mobile and thus is 
a good candidate for implementing quantum communi- 
cation schemes which are based on EPR pairs |10|. An 
important feature of this scheme, moreover, is that we 
have control over the non-locality of the entangled state: 
electron 1 and 2 are localized in different quantum dots 
and while being spatially separated from each other their 
total spin ground state (a singlet) is entangled [fljj . 

The second problem which then immediately arises 
is: How can we probe entanglement? In the following 
we provide an answer to this question and show that 
this property can be tested via mesoscopic transport and 
noise measurements. The scheme we propose consists of 
two coupled quantum dots — double-dot (DD) for short — 
which themselves are weakly coupled to two leads 1 and 
2 (see Fig. 1). In contrast to earlier set-ups involving a 
DD |I^,|l3| , we propose a scheme where an electron com- 
ing from lead 1 (2) has the option to tunnel into both dot 
1 and dot 2. This results in a closed loop, and applying 
a magnetic field, an Aharonov-Bohm (AB) phase <p will 
be accumulated by an electron traversing the DD. In the 
Coulomb blockade (CB) regime we find that due to 
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FIG. 1. a) Double-dot (DD) system containing two elec- 
trons and being weakly coupled to metallic leads 1,...,4, each 
of which is at the chemical potential The tunneling 

amplitudes between dots and leads are denoted by F, T. The 
tunneling (t) between the dots results in a singlet ground 
state with energy E s and a triplet state with Et, with 
J = Et — E s describing the effective exchange coupling be- 
tween the electron spins (qubits). The closed tunneling path 
between dots and leads 1 and 2 encloses the area A. b) Energy 
diagram of the DD system described in a) in the CB regime 
where cotunneling is dominant, i.e., U > ± > J, with 
U being the Coulomb repulsion within a single dot. 

cotunneling fl4|| the current and the noise depend on 
the state of the DD: the AB oscillations for singlet and 
triplets have opposite sign. The amplitude of the AB 
oscillations provides a measure of the phase coherence 
of the entangled state, while the period -via the en- 
closed area-provides a measure of the non-locality of the 
EPR pairs. The triplets themselves can be further distin- 
guished by applying a directionally inhomogeneous mag- 
netic field which adds a Berry phase leading to beating. 
Finally, we also discuss finite frequency noise, and show 
in particular that the power spectrum for a DD contains 
a non-vanshing odd-frequency part which is sensitive to 
entanglement. 

Model System. The DD system (see Fig. 1) contains 
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4 metallic leads which are in equilibrium with associated 
reservoirs kept at the chemical potentials fa, i = 1, . . . , 4, 
where the currents U can be measured. The leads are 
weakly coupled to the dots with amplitudes T and T, and 
the leads 1, 2 are coupled to both dots and play the role of 
probes where the current and noise are measured. The 
leads 3 and 4 are feeding electrodes to manipulate the 
electron filling in the dots. The quantum dots contain one 
(excess) electron each, and are coupled to each other by 
the tunneling amplitude i, which leads to a level splitting 
HQ J = E t — E s ~ At 2 /U of the single-particle energy 
levels in the DD, with U being the single-dot Coulomb 
repulsion energy, and E s / t are the singlet/triplet energies 
p5| . We recall that for two electrons in the DD (and 
for weak magnetic fields) the ground state is given by a 
spin singlet ||,[^] . For convenience we count the chemical 
potentials fa from E s . The coupling T to the feeding 
leads can be be switched off while probing the DD with 
a current. From now on we assume that T = unless 
stated otherwise. 

Using a standard tunneling Hamiltonian approach [fl6f , 
we write H = H + V, where the first term in H = Hfj + 
Hi + H 2 describes the DD and the leads (assumed 
to be Fermi liquids). The tunneling between leads and 
dots is described by the perturbation V = V1 + V2, where 

V n = T^2[Dl a c n ,, + ti ta D n ,.], 

S 

D n , a = e ±tlfi / 4 d hs + e^/ 4 d 2tSl n = l,2, (1) 

and where the operators c ntS and d n>s annihilate electrons 
with spin s in the nth lead and in the nth dot, resp. The 
Peierls phase tp in the hopping amplitude accounts for 
an AB or Berry phase (see below) in the presence of a 
magnetic field. The upper sign belongs to lead 1 and the 
lower to lead 2. Finally we assume that spin is conserved 
in the tunneling process. For the outgoing currents we 
have I n — ieT^2 s [D n s c n , s — c n s Z?„ iS ] . The observables 
of interest are the average current through the DD sys- 
tem, / = {I2), and the symmetrized cross correlations 
(noise) of the outgoing currents, S(t) = Re(Sl2(t)SIi(0)) , 
where SI2 = I2 — I etc. First, we evaluate the transport 
current and then the noise. 

Cotunneling current. From now on we concentrate on 
the CB regime where we can neglect double (or more) 
occupancy in each dot for all transitions including virtual 
ones, i.e. we require fa.2 < U. Further we assume that 
fii 2 > J, ksT to avoid resonances which might change 
the DD state (see also below). The lead-dot coupling T 
is assumed to be weak so that the state of the DD is not 
perturbed; this will allow us to retain only the first non- 
vanishing contribution in T to / and S(t). Formally, we 
require J > 2irv t T 2 , where v t is the tunneling density of 
states of the leads (see below). 

In analogy to the single-dot case we refer to 

above CB regime as cotunneling regime. In this regime 
electrons tunnel one by one through virtual states of 
the DD (see Fig.l). For real quantum dots not all of 



our assumptions might be perfectly satisfied and phase 
coherence might be supressed. Still, since it will turn 
out that in the cotunneling regime the current has a 
phase-coherent part with AB oscillations, it should be 
possible to extract this part even if its amplitude is 
much smaller than the incoherent one (this is a com- 
mon situation in mesoscopic transport experiments) |l9| . 
Thus, from now on we will concentrate on cotunneling 
only. We have specified now all assumptions |T^] un- 
der which the following results are valid, again they are 
U > /ii, 2 > k B T,J> 27w t r 2 . 

Continuing with our derivation of /, we note that the 
average (...) = Trp{. ..} is taken with respect to the 
equilibrium state of the entire system set up in the distant 
past before V is switched on |flq] . Then, in the interaction 
picture, the current is given by 



I = (U^I 2 (t)U), U = Texp 



dt'V(t') 



(2) 



The leading contribution in T to the cotunneling current 
involves the tunneling of one electron from the DD to, 
say, lead 2 and of a second electron from lead 1 to the 
DD (see Fig.l). This contribution is of order V 2 V 2 , and 
thus / cx r 4 , as is typical for cotunneling jl4j . 

After some manipulations we can express the current 
(||) in terms of spinless Green functions of the leads eval- 
uated at coinciding points, Gf(t) = — i(cj(t)cj(0)) and 
Gf{t) = i(4(0) Cl (t)) @|. The expanded expression for 
/ is too lengthy to be presented here, so we only mention 
that all terms of the form G± G 2 contain divergencies at 
small energy transfer around the DD levels due to reso- 
nances |Q . This is seen from the structure of Gf : first 
an electron is created in the leads and then annihilated, 
thus describing a resonant current away from the DD 
into leads 1 and 2. If the state in the DD is maintained 
through the coupling to feeding leads 3 and 4 Q , these 
resonances will be smeared and the divergencies have to 
be cut at an energy e ~ 27rz/ t r 2 . On the other hand, 
in the cotunneling regime defined above we are far away 
from such resonances, and such a resonant current is sup- 
pressed, and all divergent terms of the form G^G 2 van- 
ish. We arrive then at the following compact expression 
for the cotunneling current 

/ = C E Pi\{i\D{ s ,D 1 , s \f)\ 2 [Fi2-F2i], (3) 



Fia = J ^G<(e + E f )G>(e + E t 



(4) 



For simplicity, we assumed that the upper and lower 
branches of the tunneling loop (Fig.l) are identical. Eqs. 
(||) and ([|) show that in the cotunneling regime the ini- 
tial state \i) (with weight pi and energy Ei) of the DD is 
changed into a final state |/) (energy Ef) by the travers- 
ing electron. However, due to the weak coupling T, the 
DD will have returned to its equilibrium state before 
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the next electron passes through it [Q. For small bias, 
I Mi — 1 < J, only elastic cotunneling is allowed [fj2[ , 
i.e. Ei = i?/. However, this regime is not of interest 
here since singlet and triplet contributions turn out to 
be identical and thus indistinguishable |23|. We thus fo- 
cus on the opposite regime, — /^l > J, where inelastic 
cotunneling |24|] occurs with singlet and triplet contribu- 
tions being different. In this regime we can neglect the 
dynamics generated by J compared to the one generated 
by the bias (" slow spins" ) , and drop the energies Ei and 
Ef in the Eq. (|J). Finally, using 1 = J2f we obtain 



eTnsfT 4 C(tp) 



Mi ~ M2 
M1M2 



C (<P) = E [( d is> d isd\ s d ls i) + cosip(d\ s ,d ls dl s d 2s i 



(5) 
(6) 



where we assumed Fermi liquid leads, G>(e) = 
-i-KV t f F {Hn - e), and G<(e) = inu t f F (e - fjt n ), where 
the tunneling density of states in the leads vt = 
— f ^ m Sp G R (e,p) is expressed in terms of the retarded 
Green function |l6| , and fp is the Fermi function. 

Eq. (|^) shows that the cotunneling current depends on 
the properties of the equilibrium state of the DD through 
the coherence factor C(ip) given in (||). The first term 
in G is the contribution from the topologically trivial 
tunneling path (phase-incoherent part) which runs from 
lead 1 through, say, dot 1 to lead 2 and back. The second 
term (phase-coherent part) in G results from an exchange 
process of electron 1 with electron 2 via the leads 1 and 
2 such that a closed loop is formed enclosing an area A 
(see Fig. 1). Note that for singlet and triplets the initial 
and final spin states are the same after such an exchange 
process. Thus, in the presence of a magnetic field B, an 
AB phase factor ip — ABe/h is acquired. 

Next, we evaluate C(<p) explicitly in the singlet-triplet 
basis: \S) = (| U) - I IT)), |T ) = ^ (| U> + | IT)), 
\T+) = I TT), and |T_) = | U). Note that only the 
singlet 1 5) and the triplet \Tq) are entangled EPR pairs 
while the remaining triplets are not (they factorizc). As- 
suming that the DD is in one of these states we obtain 
the important result 



C(<p) = 2 =p cos ip 



(7) 



Thus, we see that the singlet (upper sign) and the triplets 
(lower sign) contribute with opposite sign to the phase- 
coherent part of the current. One has to distinguish, 
however, carefully the entangled from the non-entangled 
states. The phase-coherent part of the entangled states 
is a genuine two-particle effect, while the one of the prod- 
uct states cannot be distinguished from a phase-coherent 
single-particle effect JTJJ]. Indeed, this follows from the 
observation that the phase-coherent part in G factorizes 
for the product states T± while it does not so for S,T . 
Also, for states such as | Tl) the coherent part of G van- 
ishes, showing that two different (and fixed) spin states 



cannot lead to a phase-coherent contribution since we 
know which electron goes which part of the loop. 

Finally we note that due to the AB phase the role of 
the singlet and triplets can be interchanged which is to 
say that we can continually transmutate the statistics of 
the entangled pairs S, Tq from fermionic to bosonic (like 
in anyons): the symmetric orbital wave function of the 
singlet S goes into an antisymmetric one at half a flux 
quantum, and vice versa for the triplet To. 

Next, we allow for equal population of the singlet and 
triplet states in the DD, pi = 1/4, i = 1, . ..,4. Experi- 
mentally, this can be achieved e.g. by raising the tunnel 
barrier between the dots (or by increasing the magnetic 
field) jf^]] such that J vanishes or just becomes smaller 
than IcbT. Tracing then over the 4 terms we see that the 
EPR pairs, S, To, cancel each other while the unentangled 
terms, T±, add up. [This tracing can also be performed 
over the standard product spin basis as typically used for 
Fermi liquids.] Thus, in this case we return to the usual 
Fermi sea situation where we can no longer distinguish 
single- from two-particle phase-coherence j2j|. We note 
that this effect can be exploited to search for entangle- 
ment: first prepare the DD in its entangled ground state 
(singlet) with \fj,\ ±M2| > J > ^bT, and then reduce J 
below ksT; after some spin relaxation time an equal pop- 
ulation of singlet and triplet states (all contributing with 
weight 1/4) is reached with a concommitant increase of 
the current by a factor of 5/2. 

We would like to stress that the amplitude of the AB 
oscillations is a direct measure of the phase coherence of 
the entanglement, while the period via the enclosed area 
A = h/eBo gives a direct measure of the non-locality of 
the EPR pairs, with Bq being the field at which <p = 1. 
The triplets themselves can be further distinguished by 
applying a directionally inhomogeneous magnetic field 
(around the loop) producing a Berry phase <& s p6| , 
which is positive (negative) for the triplet m = 1(— 1), 
while it vanishes for the EPR pairs S, To. Thus, for J = 
we will eventually see beating in the AB oscillations due 
to the positive (negative) shift of the AB phase <f> by the 
Berry phase, ip = <&±$ B . We finally note that the closed 
loop shown in Fig.l can actually be made as large as the 
dephasing length by replacing the dotted lines outside 
the dots e.g. by wave guides forming a loop with a lead 
attached somewhere to it f23| . 

Shot noise. We evaluate now the cross correlations S(t) 
of the currents by expanding in powers of V, see (||). The 
first non- vanishing contribution is of order V\V-z. Under 
the same assumptions as before, we calculate the spectral 
density of the noise, S(u) = J dte lu]t S{t). For the zero- 
frequency noise we obtain 



S(0) 



e 2 r 4 



2tt 



E 

hf,s,s' 



Pi \(i\Dl.D 1<s \f)\ 2 [F 12 +F 21 ]. (8) 



We compare now to Eq. (||) and note that depending on 
the sign of \i\ — p.2 either F\i or F 2 \ vanish (at T = 0). 
Therefore, in the cotunneling regime the noise assumes 
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its Poissonian value, 5(0) = — e|J|, and we see that the 
Fano factor (noise-to-current ratio) is universal and the 
current and its cross-correlations contain the same in- 
formation. [When we allow for resonances of the type 
discussed above, this is no longer the case [^3) .] 

For finite frequencies in the regime \fii — fi2\ > J we 
obtain after lengthy calculations, 

S(u) = (eirvtV 2 ) 2 [X u + Xl u ] , 
C(<p) 
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The real part of 5(o>) is even in ui, while the imaginary 
part is odd. A remarkable feature here is that the noise 
acquires an imaginary (i.e. odd-frequency) part for finite 
frequencies, in contrast to single-barrier junctions, where 
Im5(a->) always vanishes since we have 5I± = —SI2 for 
all times. In double-barrier junctions considered here we 
find that at small enough bias A/i = fii — /i 2 <C \i = 
(fi± + H2) /2, the odd part, Im5(w), given in (|9|) exhibits 
two narrow peaks at u) = ±/i, which in real time lead to 
slowly decaying oscillations, 



Sodd{t) = e 2 ™ t 2 r 4 c(^) 



sin(A/it/2) 
fit 



sin(/xt). (11) 



These oscillations again depend on the phase-coherence 
factor C with the same properties as discussed before. 
These oscillations can be interpreted as a temporary 
build-up of a charge-imbalance on the DD during an un- 
certainty time ~ juT , which results from cotunneling of 
electrons and an associated time delay between out- and 
ingoing currents. 

Finally, allowing currents to flow from feeding leads 
via DD into leads 1, 2 we can arrange for scattering of 
unentangled electrons (as considered previously in noise 
p7[) but now also of entangled ones. In the latter case 
we get a non-trivial Fano factor [|8| due to antibunching 
(triplets) and bunching (singlet) effects in the noise p0| . 
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